function w = Aop(v)
    w = v;
    w(2:end-1) = 2 * v(2:end-1) - v(1:end-2) - v(3:end);
end